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ABSTRACT 


Three numerical model simulations of cyclogenesis are compared to examine the role 
of boundary laver stratification in enhancing the components forcing the ageostrophic 
circulation of an idealized warm front. The Sawver-Eliassen diagnostic equation 1s ap- 
plied to examune the contributions of frictional forcing and diabatic heating, as Well as 
confluent and shear geostrophic deformation, to forcing of the secondary circulation of 
a warin front. The surface heat and moisture flux distributions are varied in each case, 
in order to evaluate the effect on cach component as well as the cyclogenesis. 

Results confirm previous studies that geostrophic deformation forces strong 
frontogenesis at the surface, and at mid-levels frontogenesis is Weaker and forced pri- 
marily by latent heat release. Frontogenctical forcing is modified by small-scale fric- 
tional forcing and diabatic heating, which depend upon the surface and boundary laver 
processes. Although frictional forcing comprises less than five percent of total forcing 
of the Warm front. it apparently enhances frontogenesis, partially due to indirect effects 
on other more dominant processes. The intensity of frictional forcing is strongly de- 
pendent on the surface heat flux distribution and track of the evclone relative to the 
sea-surfuce temperature gradient. In the absence of surface fluxes, frictional forcing 1s 
negligible. Surface forcing due to diabatic heating is frontolytical. and reduces total 
forcing at low levels by about twenty percent. The magnitude of the frontolstical forcing 
by digbatic heating at low levels is only partially dependent on the surface heat and 
mieisiure flux distribution. Even in the absence of surface heat and moisture fluxes, the 
frontolytical forcing persisted, suggesting problems with the cumulus parameterization 
or errors in the computation of this term. These results remain to be verified against an 


actual case study of a warm front. 
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I. INTRODUCTION 


A. CYCLOGENESIS AND FRONTOGENESIS 

Mid-latitude cyclones are transicnt waves in the global zonal flow that perform an 
important role in transporting heat and momentum in the atmosphere. Historically, 
cyclone development was explained as a wave instability on the polar front. with the 
existence of the front considered a necessary precursor for cyclogenesis. This theory, 
known as the polar front theory of cvclogenesis, significantly influenced Petterssen 
(1956) who related cvclogenesis to the interaction of vorticity advection aloft with sur- 
face frontal features. The modern view, based on baroclinic instability theory, considers 
fronts to be a consequence of the development of a cyclone and a source of baroclinic 
energy for cyclogenesis. The baroclinicity of the atmosphere and the synoptic-scale flow 
combine to lead to cvclogenesis, the development of sub-synoptic scale fronts, and the 
formation of jet streaks. 

‘The basic prenuse of baroclinic instability theorv is that the circulation of a cyclone 
re ults from any small perturbation forced by advection of positive (cyclonic) vorticity 
aloft with the advance of the upper-level wave or by warm advection at low levels over 
a broad region. Gcostrophic deformation associated with cvclogenesis acts to enhance 
the temperature gradient in particular regions of a cyclone. This process is known as 
frontogenesis. Vertical and divergent ageostrophic motions arise in response to 
lrontogenesis, and are the mechanism for the adjustment of vorticity and thermal prop- 
erties to return the svstem to thermal wind balance. The agcostrophic circulations as- 
sociated with frontogenesis are small-scale manifestations of the basic synoptic-scale 
forcing. Although these synoptic-scale processes are thought to be the primary proc- 
esses for cvclogenesis, frontogenesis is a secondary process that may feed back on the 
synoptic-scale forcing to produce stronger cyclogenesis than that given by just the 
synoptc-scale forcing. The interaction between svnoptic-scale vorticity and temperature 
advection and mesoscale frontogenesis in cyclones is not vet completely understood. 

I. Sub-Svnoptic Scale Processes 

Recent cyclogenesis studies suggest important coupling between cyclogenesis 
and mesoscale processes associated with frontogenesis and jet streaks. Ucccllini et al. 
(1985, 1987) and Keyser and Shapiro (1986) have shown that small scale upper-level 


baroclinic processes associated with upper-level jet streaks enhance the vertical circu- 


lation and cvclogenesis. Strong jet streak divergence induces a secondary circulation 
that extends through the depth of the troposphere. This enhances the development of 
a low-level jet and cvclogenesis at the surface. This mesoscale process is considered to 
be an important mechanism for some explosive cyclones. 

Another important mechanism in some explosive cyclones is the release of latent 
heat in the mid-troposphere. The release of latent heat also influences the frontogenesis 
process. Wilhams etaal. (1981), using a numerical model, and Thorpe (1984). using an 
analytical quasi-geostrophic model, showed that latent heat release on the warm side of 
the front is frontogenetical aloft. The heating directly increases the temperature gradient 
at levels of large heating. and the induced divergent circulation leads to convergence 
below the heating which increases the temperature gradient in that region. The precipi- 
tation maximum occurs where the vertical motion is a maximum. Evaporative cooling 
occurs When precipitation falls through drv air below on the cold side of fronts, thus 
enhancing the cross-front temperature gradient ahead of the front. 

2. Planetary Boundary Layer Processes 

Although baroclinic synoptic and mesoscale processes mentioned above are 
dominant, the importance of surface and boundary laver processes also has been dem- 
onstrated in recent studies. Fhe boundary laver processes contribute to cvclogenesis and 
frontogenesis through thermodynamic effects due to the heat and moisture fluxes, and 
dvnanne effects due to friction. [Friction further disrupts geostrophic balance and in- 
creases the ageostrophic component of the wind through frictional wind turning. This 
results in an additioual component of convergence and shear due to friction. Kevser and 
Anthes (1982) showed that convergence due to surface friction acts frontogenetically at 
the surface. Behind the frontal zone the vertical shear of the ageostrophic wind, which 
is frictionally induced, acts to diminish the static stabilitv in the planetarv boundary 
laver (PBL) and enhance the static stability at the top of the boundary laver. Ahead of 
the front. differential temperature advection by the agcostrophic component of the wind 
acts to increase the static stability in the boundary laver and decrease the static stability 
at the top of the boundary laver. The result is a compressed thermal gradient in the 
PBL. which is frontogenetical. Keyser and Anthes (1982) also showed that frictional 
convergence drives a vertical velocity jet above the suriuce fromt. This vertical jet is 2 
component of the secondary circulation. 

Danard and Ellenton (1980) evaluated the surface heating in several east coast 
U. S. evclones, and concluded that the surface heating acted to decrease the low-level 


baroclinicity by warming the cold advection regions and vice versa. This finding is 


tl 


consistent with previous studies of the effects of surface heating in baroclinic instability 
models, and indicates a frontolvtical role of the surface heat flux. However, in studies 
of the Presidents’ Dav storm of 1979, Bosart (1981) and Bosart and Lin (1984) indicated 
that the low-level contributions by the surface heat and moisture fluxes were important 
to cvclogenesis. They found that positive heat fluxes to the northeast of the surface low 
potentially contributed to rapid development through distribution of the cyclogenctic 
environment. Reed and Albright (1986) found similar results in their study of an eastern 
North Pacific cyclone. 

\ favorable surface heat and moisture flux distribution can also enhance 
frontogenests and ecvclogenesis through the associated latent heat release above the 
boundary layer. Chen et al. (1983) found that mid-tropospheric latent heating decreased 
when the surface heat and moisture fluxes were removed. Atlas (1987) found that the 
vertical transport of air that is heated and moistened bv the surface heat and moisture 
fluxes can result in significant mid-tropospheric heating through the release of latent 


he ut. 


fae ix e CYCLONES 

Compared with continental and coastal cyclones, marine cyclones have a more 
complex surface heat and moisture flux distribution contributing to frontogenesis and 
evclogenesis, due to the influence of the sca-surface temperature (SST) distribution. 
Businger and Shaw (1984) suggested that for a simple uniform flow across an SST gra- 
dient, PBL convergence will result from stratification effects on the frictional wind 
turning as the wind flows across an SST gradient. As a evclone moves across a fixed 
SST distribution, changes in the frictional convergence and the transport of boundary 
laver air to the mid-troposphere are likely to occur. In their study of surface fluxes and 
divergence in Storm Transfer and Response Experiment (STREN) storms, Fleagle and 
Nuss (1985) showed that the distributions of the frictional convergence and vertical ve- 
locity at the top of the PBL are primarily determined by the horizontal distribution of 
the surface wind stress. In turn, the wind stress depends not only on the synoptic-scale 
forcing at the top of the PBL, but also on the PBL stratification, the surface heat flux 
distribution, and thermal wind turning in the PBL. Nuss (1989) showed that stable 
stratification behind the warm front, and unstable stratification ahead of the warin front 
in a region of upward heat fluxes, apparently resulted in enhanced frictional convergence 


ahead of the front. Danard (1986) found that increasing the SST by five degrees resulted 


in greater boundary laver convergence and enhanced convective heating in the model 
cyclogenesis. 

Nuss and Anthes (1987) studied the effeet of surface fluxes and boundary laver in- 
fluences on an idealized marine cyclone. They showed that although baroclinic insta- 
tility, which depends on the vertical wind shear and static stability, as shown originally 
by Charney (1947) and Eady (1949), is sufficient for rapid cyclogenesis, relatively small 
changes in the low-level static stability, meridional temperature gradient, and surface 
heat and moisture fluxes significantly modified the development of the model cyclone. 
These low-level changes resulted entirely from varving the SST distribution. They found 
that the deepening rate of the cyclone increased by 25% for complete removal of the 
surface heat and moisture fluxes, and increased bv 15% for a modified SST distribution 
in phase with the low-level temperature pattern. 

These results were further anaivzed by Nuss (1989), in order to clarify the role of 
surface fluxes in the overall development of the cvclone. He found that the effect of 
surface heat and moisture fluxes on the structure and dyvnamies of oceanic cvclones can 
be substantial and that strictly boundary layer processes act to enhance the vertical cir- 
culation and evclogenesis through a dynamic feedback process. Specifically. he found 
that the effect of surface heat and moisture fluxes depends critically on their distribution 
relative to the cvelone. For a zonal SST distribution. patterns of surface heating 
counteract BL temperature advection and oppose cyclogenesis, as found in previous 
studies. including Danard and Ellenton (1980). For a sinusoidal SST distribution es- 
senuallv in phase with the low-level temperature distribution, the resultant pattern of 
surface heating vielded unstuble stratification to the northeast and stable stratification 
south of the warm front. This acted to enhance cyclogenesis through enhanced frictional 
convergence as well as increased heat and moisture transport to the surface cyclone from 
the northeast. 

The relationship of enhanced frictional convergenee to other frontogenesis processes 
Was not examined and is the focus of this study. Although the surface heating pattern 
may be a damping effect, associated mid-tropospheric latent heat release increases the 
temperature gradient and cyclogenesis, which indireetly enhances the PBL frontogenesis 
processes. This indirect relationship between surface diabatic processes and the intensity 
of the geostrophic frontogenesis and cyclogenesis is not well understood. Uccellini et 
al. (1987) have demonstrated the importance of the non-linear coupling between 
boundary layer processes, latent heat release. and jet stream forcing along a coastal front 


for the Presidents’ Dav storm. For open ocean cyclones, the warm front region 1s likelv 


to be the focal point of the interaction between surface fluxes, diabatic heating, and the 
jet, as suggested by Nuss (1989), and will be the subject of this study. Studies by 
Emanuel (1985, 1988) and Sanders (1986) also suggest that latent heat processes associ- 


ated with the warm front may contribute to enhanced cyclogenesis over the ocean. 


Cy OBJECTIVE 

The purpose of this research is to quantify the contributions of various forcing 
mechanisms to the secondary circulation around a warm front, and assess their impact 
on the development of a cyclone. The influence of the surface heat and moisture fluxes 
on the frictional forcing and diabatic heating, and resultant secondary circulation, is of 
particular interest since it is not as clearly understood as the larger-scale geostrophic 
deformation processes. In addition, the feedback of these diabatic processes on the 
geostrophic deformation and warm frontal intensification will be examined. The specific 


Gopeetivcs arc: 


1. Diagnose the geostrophic versus diabatic processes in warm frontogenesis 


to distinguish frictional convergence from pure geostrophic deformation; 


tu 


Describe the feedback mechanism of surface heat fluxes on geostrophic and 
Jiabatic forcing by quantifving their frontogenetic effects on the friction 
and latent heat release and their frontolvtical effects on the PBL thermal 


gradient: 


3. Describe the time evolution of the forcing: 


4. Obtain an understanding of the interaction of the SST with fronts, 


The Sawver-Eliassen diagnostic equation will be used to diagnose the forcing of the 
vertical circulation around the warm front in three numerical simulations of 
evclogenesis. The Sawver-Eliassen equation is reviewed in Section II, followed bv a brief 
description of the model experiments in Section III. Section IV compares the results 
of the three simulations. A discussion of these results is given in Section Y, followed by 


conclusions and recommendations for future research in Sections VI and VII. 


Il. BACKGROUND 


Sawyer (1956) showed how the secondary circulation around fronts can be calcu- 
lated on the basis of quasi-geostrophic theory, and Eliassen (1959, 1962) extended his 
work to include friction and latent heating effects, as well as the vertical shear in the 
cross-front plane, to the forcing of the secondary circulation. They developed the theory 
of describing temporal changes in the vertical and horizontal gradients of temperature 
and absolute momentum in the cross-front plane in terms of the induced geostrophic and 
ageostrophic deformation field. It 1s from this theory that the Sawyer-Eliassen diag- 
nostic equation for the secondary circulation was derived, This equation enables diag- 
nosis of small-scale diabatic and frictional processes as Well as those processes associated 
with geostrophic deformation forcing the secondary circulation around a front. Ehassen 
(1962) demonstrated that the circulation is counter-clockwise (cyclonic) around regions 
of positive forcing, and clockwise (anti-cvclonic) around regions of negative forcing. The 
circulations are concentric ellipses mn the vicinity of the point source of forcing (Fig. 1). 
oriented verucally along the front and becoming vertically or horizontally distorted at 
some distance from the pomt source as a function of vorticity and mertial stability 
(Kevser and Shapiro 1986). The Sawver-Llhassen equation was later expanded by 
Hoskins and Draghici (1977) and Hoskins et al. (1978) to account for variations m the 
along-front plane, thus generalizing the theory of secondary circulation into three di- 
mensions, known as the Q-vector theory. The resultant eyuation is known as the 


three-dimensonal (3-D) form of the Sawver-Lhiassen equation. 


A. AGEOSTROPHIC CIRCULATION AROUND A WARM FRONT 

Gidel (1978) describes the processes of geostrophic balance and continuity of mass 
flow which lead to the secondary circulation around a warm front. The rising motion 
within the front coupled with the sinking motion on the cold air side of the front is a 
thermodvnamically direct circulation. By continuity of mass, the flow ts convergent at 
low levels in the upward vertical motion region and divergent at upper levels. Transverse 
flow results at lower and upper levels, completing the circulation cell in the region. 
Upper-level westerly (positive) momentum and low level easterly (negative) momentum 
are created, which contribute to increasing the vertical shear of the along-front wind 
component. The adiabatic cooling of the rising warm air and adiabatic warming of the 


sinking cold air contribute to decreasing the temperature gradient across the front. Figs. 





Fig. 1. Streamlines around a positive point source. Streamiimes in the (a) vy-plane 


and in the (b) sp-plane. 21 denotes absolute momentum. From Ehassen 
(1962). 


2 and 3, from Gidel (1978), depict the initial and advected potential temperature surfaces 
and 1sotachs, and the induced secondary circulation around a warm front. In the frontal 
zone the wind field is stronger than that required for geostrophic balance, implying that 
Coriolis dominates the PGT, which sets up a transverse flow to the right (north). Im- 
mediately ahead of the front, the opposite is true, and the resultant transverse flow is to 
the left (south). The flow is convergent directly beneath the front. and divergent behind 
and ahead of the front. The convergence near the leading edge of the frontal zone acts 
as a deformation field which enhances the potential temperature gradient across the 
front. The circulation cell behind the front is thermodynamically indirect and is the in- 
direct cell discussed by Eliassen (1962) with respect to warm fronts. In the indirect cir- 
culation, the cross-front thermal gradient is enhanced by adiabatic processes and the 
along-front vertical wind shear is diminished by low-level divergence and upp .-r-level 


convergence. 
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Idealized cross section through a warm front. The initial position of the 


potential temperature surfaces is indicated by the solid lines, and their po- 
sition after they have been slightly advected is indicated by dotted lines. 
The dashed lines are isotachs of the along-front wind field in geostrophic 


balance with the initial potential temperature field. Irom Gidel (1978). 
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Fig. 3. Secondary circulation around a warm front. The dotted lines are repeated 
from Tig. 2. The dashed lines are the isotachs of the along-front wind field 
after being advected from their initial position shown in lig. 2. The solid 
lines are the isotachs of the along-front wind field « which geostrophically 
balances the potential temperature field, shown as dotted lines. The heavy 
arrows indicate horizontal and vertical components of the induced cireu- 
lation. From Gidel (1978). 


The ageostrophic flow around a positive center of forcing ts counter-clockwise, 
which gives the sense of the secondary direct circulation cell if it was located ahead of 
the warm front. Similarly, a negative (clockwise) center gives the direction of flow of the 
indirect cell if the center was located behind the front. Fig. 4a is a depiction of the the- 
oretical location of the positive and negative ellipses of forcing around a warm front. 
A dipole effect is set up between the two centers, and the gradient of forcing from one 
cell to the other is depicted as strongest across the vertical flow regions, assunung the 


downward and upward motions are the same magnitude. Theoretically, the stronger the 


magnitude and gradient of the cells, the stronger the magnitude of flow of the circulation 
cell. If the updraft region ts stronger than the downdraft region, one would expect the 
gradients to be as in lig. 4b, 

The above discussion assumes the scale of geostrophic and diabatic forcing is the 
same as that of the circulation cell. However, processes which contribute to enhance- 
ment of the circulation cell are not on the same scale. The Sawyer-Eliassen equation is 
a means by which the small-scale processes of frictton and diabatic heating can be in- 
cluded with geostrophic deformation to diagnose the forcing of the cell. Total forcing 
would not necessarily be a single svmmetric center, but rather a center distorted by 
smaller regions of positive and negative forcing acting to enhance or retard the second- 


ary circulation. For example. one would expect enhanced regions of forcing at low levels 


duc to geostrophic and frictional confluence, and at mid-levels duc to latent heat release. 





Fig. 4. Ageostrophic forcing centers around a warm front. Lllipses of forcing are 
for (a) an updraft and downdraft of equal strength, and (b) a stronger up- 
draft. Note the ortentation of the front and x- and y-axes have been re- 


versed from Figs. 2 and 3. 


B. SAWYER-ELIASSEN EQUATION 
Lhassen (1962) introduced the quantity absolute momentum to describe the second- 
ary circulation in terms of an absolute frame of reference, the equations for which can 


be written as 


Mm. =u, — fr (1) 
and 
i ee ey (2 


The thermal wind relationships in pressure coordinates are 
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Following Hoskins and Draghici (1977), and Kevser and Shapiro (1986), prognostic 


equations for absolute momentum and thermal advection can be written as 
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By rearranging the terms, the x and y components of the Sawver-Ehassen equation 


can be written as 
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By taking the dot product of —A XN +V@ and the vector equation made up of the x 
and y components given above, the magnitude of the forcing and circulation can be 


obtained. The nght hand side (forcing) of the equation becomes 
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The terms are then normalized by dividing by | kN ;Vé@ |. to obtain the magnitudes 
of the components contributing to total kinematic forcing of the warm front. The 
Sawver-Lhassen equation is not solved: rather this form enables a diagnosis of the con- 
uibuuons to the forcing of the circulauon cell by each term. Centers of positive and 
negative forcing are derived by this diagnosis, and represent centers of counter-clockwise 
and clockwise circulation respectively. The 3-D form is used in place of the 2-D version 
because the data set is three-dimensional and because the along-front gradients are not 


assumed to be in geostrophic balance. 


€C. DESCRIPTION OF TERMS 
1. Geostrophie Confluence 

The first term in equation (11) represents the effect that confluent geostrophic 
deformation has in enhancing the horizontal temperature gradient (V7) by compressing 
the isotherms across the front. The strongest centers of confluence (positive forcing) and 
diffluence shouid be as in Fig. $8. The low-level confluence and upper-level diffluence 
represent a downward transport of momentum, which acts to reduce the vertical wind 
shear. while at the same time the thermal gradient is increasing. Following Kevser and 

ef 


P : : ae ae, 
Shapiro (1986), the resultant contributions to Nee) / dt and dl; 5, ) / dt are respec- 
Al: 


tively negative and positive, Which is a tendency towards imbalance in the thermal wind. 





a frontogenetic process. This imbalance is counteracted by the secondary direct circu- 


lation :orced by the confluence, which simultaneously increases the verucal wind shear 


Me 


through the net upward transport of positive momentum, while reducing the temper- 
ature gradient by adiabatic processes. 
2. Geostrophic Shear 

The second term represents the effect that gcostrophic shear deformation has in 
rotating the along-front temperature gradient perpendicular to the front. If the winds 
across the front increase to a maximum at the low Icvel easterly jet on the cold air side 
(ahead) of a warm front, cvclonic shear is present. I{ the easterly wind ahead of the front 
is Stronger than that behind the front, differential cyclonic shear in the presence of 
along-front decreases in potential temperature tends to enhance the cross-front temper- 
ature gradient. Larger negative values of u, ahcad of the front relative to weaker values 
behind the front (st <0) diflerentially advect higher values of potential temperature 
at a greater rate north of the front relative to south of the front. At the same time, the 
oY 
cp 
shear is cyclonic,’ thus reducing the magnitude of the vertical wind shear, creating an 





vertical wind shear (=—) tends to tip the , isotachs into the vertical where the horizontal 
imbalance. Again, the direct circulation would counteract the imbalance created bv the 
shear as with the confluence. One would expect that the strongest forcing due to shear 


would be between the low-level casterly jet and the front. 
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Fig.=3: Transverse ageostrophic circulations for an idealized front. Circulations 
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oa 0, a > 0). Dashed lines, isotachs of w, (denoted by U); thin solid 
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lines, isotachs of y, (denoted by VY); thick solid lines, streamlines of 

ageostrophic circulation. Adapted from Lhassen (1962). From Keyser and 


Shapiro (1986). 


3. Friction 

The third term represents the convergent effect that friction has on the 
geostrophic wind. Winds at the surface experience friction, but are maintained in part 
by the downward momentum flux from stronger winds above the surface. Differential 
downward momentum fluxes act like confluence to compress the temperature gradient. 
As the cyclone intensifies the winds increase. which leads to increased friction and fric- 
tonal turing of the winds toward the low. This enhances frictional convergence. Ifthe 
winds were to increase uniformly around the low, there would be no differential mo- 
mentum flux across the front. However, the wind gradient across a frontal zone 1s 
greater than elswhere around the low, so the maximum differential friction is in the vi- 
cinity of the fronts. 

As discussed in Section I, favorable surface heat fluxes also enhance frictional 
convergence. Regions of upward heat flux into the PBL from a relatively warm surface 
can be likened to warm air advection into the region. This would cause a vecring of the 
PBL wind with height, or an increased turning angle and further convergence at the 
surface. Similarly, a downward heat flux from the PBL to a cool surface would cause a 
backing of the thermal wind with height, or a reduced angle of turning and reduced 
convergence. Tor a warm front in which there are downward heat fluxes in the cold air 
sector ahead of the warm front as well as in the warm air sector behind the front. one 
can expect that for a given SST distribution, the air-sea temperature difference would 
be greater in the warm air region behind the front than in the cool air region ahead of 
the front. The downward heat flux would be greater behind the front than ahead, with 
the backing of the thermal wind stronger behind the front than ahead. This creates a 
net differential across the front, with stronger convergence ahead of the front. Thus, the 
frictional convergence is composed of the geostrophic component and a heat flux com- 
ponent. Fig. 6 depicts the geostrophic frictional forcing as a result of the differential 
downward transfer of momentum from the winds above the PBL, and the heat flux 
component of frictional forcing within the PBL duc to surface heat fluxes. The depiction 
is for the case of stronger downward momentum fluxes behind the front than ahead, and 
for stronger downward heat fluxes at the surface behind the front. This would result in 
a reinforcement of frictional confluence. The greatest differential would be across the 
warm front, and the strongest frictional confluence would be located where the strongest 
differential occurs. Realistically it 1s possible to have stronger downward momentum 
fluxes ahead of the warm front, which would oppose the heat flux component depicted 


above, resulting in neghgible frictional forcing. 
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Fig. 6. Centers of frictional forcing. [leavy lines above PBL denote secondary 


circulation cells, with vertical arrows in the PBL denoting downward mo- 
mentum fluxes (solid) and downward heat fluxes (curved). Leavy arrows 
indicate stronger downward momentum fluxes. Dashed Imes depict warm 


front. 


4. Diabatic Heating 

The fourth term represents the effect diabatic heating has on the horizontal 
temperature gradient. In the boundary laver, differential leat Mixes act to relax tie 
temperature gradient, and are strongly dependent upon the surface temperature distrib- 
ution, as discussed in Section J]. Surface heat and moisture fluxes in the PBL can be 
turbulently mixed out of the PBL. providing a source of heat and moisture to the mid- 
troposphere in addition to the heat and moisture advected mto the region. At mid- 
levels, latent heat release within and behind the warm front, and evaporative cooling 
ahead of the front enhance the temperature gradient. One should expect centers of 
positive and negative diabatic forcing to be located as in Vig. 7. The low level forcing 
opposes the secondary circulation at low levels and the mud level forcing enhances the 


circulation. 





big: 7. Centers of diabatic forcing. Forcing in the PBI opposes the secondary 


circulation. In the mid-troposphcre it enhances the circulation. 


Ly 


III. DESCRIPTION OF DATA SET 


A. NUMERICAL MODEL 

The numerical simulations for the three case studies used in this study were run for 
48 hours using the Pennsylvania State University National Center for Atmospheric Re- 
search (PSU NCAR) mesoscale model version 4 (M4), which is described in detail by 
Anthes et al. (1987). The usual metcorological data plus the boundary layer fluxes were 
output everv 3 hours. The version of the model used for this study utilized 14 vertical 
lavers from the surface to 100 mb, periodic east-west boundaries, a # - plane variation 
of the Coriolis parameter, and the high resolution boundary laver model. The boundary 
laver occupies five of the 14 layers, with the lowest laver within the surface laver. 

Iiiferential friction can result from the boundarv layer structure and processes in 
the model since the model vertical heat. moisture, and momentum fluxes use a vertical 
diffusion coefficient that depends on the local bulk Richardson number (R,). The 
boundary conditions are considered stable when the bulk Richardson number (R,) is 
greater than zero and Z,// is large. Z, 1s the height of tle stirlace layer and lL is the 
Monin-Obukhovw length. Tor very stable conditioms (&, = 0.2 and 7./L unrestricted), 
10 turbulent fluxes oceur. Tor stable conditions (0.0 < R,< 0.2 and Z,/L unrestricted), 
damped meehanical turbulent fluxes occur, and for marginally unstable or neutral con- 
ditrons (R, < 0.0 and Z,/L < 1.5). both mechanical and buoyant turbulent fluxes occur. 
Unstable boundary conditions apply when both R,< 0.0 and Z,/L > 1.5, and under 
these conditions the vertical fluxes are computed using a convective plume model. In- 
cluded in the fluxes ts a downward entrainment flux at the top of the mixed laver equal 
to 20%of thesumigce flux: 

Cloud formation and precipitation are parameterized through both non-convective 
and convective processes. Stable lifting that results in relative humidities greater than 
100°. is treated as non-convective precipitation, and the latent heat is added directly to 
the thermodynamic equation. Unstable ascent or convection is parameterized using a 
modified Kuo (1974) scheme (Anthes 1977) which redistributes the moisture and latent 
heat vertically through prescribed profiles. The criteria for convection are total column 
moisture convergence greater than 0.0} g m2 s?} and the predicted cloud height greater 


than 700 mh 


Diabatic heating used in the Sawver-Eliassen diagnostic equation is computed as the 


residual of the thermodynamic equation 
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The time tendency was computed as a six-hour time difference and introduces some error 
into the diabatic heating values. However, the diabatic heating data does correspond 
with the rainfall regions, giving at least qualitative confidence in their distribution. Be- 
cause frictional forcing 1s computed as a second order vertical differential, forcing at the 


surface is represented by forcing at a level above the surface in the muddle of the PBL. 


B. INITIAL CONDITIONS 

The initial sea-level pressure (SLP). surface equivalent potential temperature, and 
observed wind field for the three cases are shown in ig. $8. The initial baroclinic insta- 
bility conditions are the same in each case. with the only differences due to the respective 
surface heat and moisture flux distributions. The initial model cyclone differs from an 
actual atmospheric cyclone most significantly in its lack of an initial vertical circulation, 
Its moisture structure being too drv compared with typical cloud patterns of actual 
evclones, and a sheht upshear tilt in the thermal wave. Al but the moisture distribution 
difference are chminated in the first six hours of each model simulation. The moisture 
distribuuon takes I$ hours to resemble the typical cloud pattern of an actual cyclone, 
Which also gives less certainty to the diabatic heating contributions to the model 
evclogenesis and frontogenesis. 

The control experiment (Case 1) uses a zonal SST distribution (Fig. 9) and includes 
surface heat and moisture fluxes throughout the stmulation. With this SST distribution 
the mean baroclinicity in the atmosphere is reinforced, and the pattern of surface heating 
introduces temperature tendencies that counteract the initial thermal advection (Nuss 
and Anthes, 1987). The initial surface flux distribution is characterized by upward heat 
and roisture fuxes west of the surface low and downward fluxes east of the surface low. 
The second experiment (Case 2) uses the same zonal SST distribution but insulates the 


atmosphere from the sea surface by specifving a zero heat and moisture flux condition 
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Fig. 8. Initial conditions. S1P in mb (solid), 0, in K (dashed). and observed wind 
vectors in ms for cach case. Length of wind vector denotes intensity, 


with the maximum length representing 50 m s-!. 


throughout the simulation. Because the only source of moisture im the model is from the 
surface moisture flux, this case hus no moisture directly transported to mid and upper 
levels as in the other two cases, and relics on the cloud formation and precipitation 
parameterization applicd to rismg warm air cooled through adiabatic processes. The 
downward momentum flux is not specifically modified from Case 1, except indirectly 
through stratification effects that result from no surface heating. With no surface fluxes, 
the total frictional convergence is reduced to that caused by dillerential downward mo- 


mentum fluxes, and there is no thermal damping eflect in the boundary layer. 





Fig. 9. Case 1 and Case 2 SSI distribution. SST in °C (solid), with the initial total 


(latent plus sensible) heat flux distribution in W #7 2, and the cyclone track 


and position at six-hour intervals from 00 h to 48 h indicated for Case I. 
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The final experiment, Case 3, uses a sinusoidal SST distribution and includes surface 
heat and moisture fluxes as in Case | (Fig. 10). The sinusoidal SST distribution in this 
case is in phase with the thermal structure of the cyclone, thus rcinforcing the low-level 
barochinicity instead of the large scale mean baroclinicity. ‘This distribution results in 
large upward total heat fluxes to the northeast of the surface low as well as to the west. 
In this case the magnitude of the heat flux is largest northeast of the surface low, en- 


hancing the frictional convergence northeast of the surface low compared with Case 1. 
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Fig. 10. As in Fig. 9, except for Case 3. 


C. DEVELOPMENT OF THE CYCLONE 

The differences in structure and dynanmucs of each cyclone throughout the simu- 
lations represent the direct and indirect effects of the surface flux distribution unique to 
each case, since the initial baroclinic instability factors are identical. The intensity of 
cyclogenesis in each case can be described in part by the sea-level pressure (SLP) 
deepening rate. Fig. If is a plot of the lowest closed contour of SLP around the center 
versus time for cach case. By the Sanders and Gyakum (1980) definition of explosive 
development (I mb /#-! for 24h or more), all three cyclones develop explosively, with the 
average deepening rate 1.46 mb /:- for Case 171-89 inb A“ lor Case 2, and [65 nin! 
for Case 3 (Nuss and Anthes 1987). 

In Cased, tiie cyclone tracks northeast towards colder water, tnitially over a region 
of increasing downward heat fluxes, and then towards the region of increasing upward 
heat fluxes in the top right quadrant of Fig. 9. As the cyclone moves towards this re- 


gion, the differential in static stability across the warm front is enhanced. Increased 





Fig. ii. SLP (mb) versus tine for each case. Note that the deepening rates are 


the same for cach case until 12 h, and the deepening rates of Cases 2 and 


3 are equivalent until 30h. 


upward heat fluxes ahead of the front reduce the static stability in this region compared 
with the comparatively small increase in static stability behind the front in the region 
of downward heat fluxes. Diflciential static stability in the PBL increases frictional 
convergence, however it 1s small relative to the thermal damping effects of the heat flux 
distribution, which acts to slow f1ontogenensis by warming the cold air sector ahead of 
the front and cooling the warm air sector behind the front. his prevents the formation 
of concentrated temperature gradients at the surface and associated strong temperature 
advection (Nuss and Anthes, 1987). 

The removal of surface heat and moisture fluxes in Case 2 results m an increased 
rate of cyclogenesis compared with Case I. In the absence of surface fluxes, diflerential 
Static stabilty across the front is reduced, so frictional convergence due to differential 
downward niomentum fluxes is no longer enhanced. Ilowever, the frontolvtical thermal 
effect of the Case | surface heat fluxes is removed, and indeed, Nuss and Anthes (1987) 


showed that the low-level thermal cold and warm front advection increased in this case 


compared with Case 1, and the temperature gradients were stronger, indicating stronger 
cold and warm fronts. 

In Case 3 the cyclone tracks northeast across the SST gradient, moving shghtly to- 
wards warmer Water until the 24 hour pomt (24 h) when the cyclone starts tracking to- 
wards colder water. This track results in an initial reduction in downward fluxes east 
of the surface low until 24 hours, when the downward fluxes begin to increase as the 
cyclone changes direction. The surface flux distribution of the first 24 hours favors a 
maintenance of the baroclinicity of the cyclone. The distribution of strong upward 
fluxes north of the warm front results in reduced low level static stability and stronger 
heating of the PBL in this region, which enhances {frictional convergence northeast of the 
surface low. The frictional convergence in this case is stronger than in Case 1, where the 
upward heat fluxes northeast of the surface low were not as strong. Another significant 
difference between the two cases is that the moisture flux to the northeast of the warm 
front in Case 3 is greater than 200 W wr? compared with near zero in the control case 
(Nuss, 1989), and there is an upward moisture flux throughout most of the warm sector 
in Case 3 compared with near zero fluxes in Case 1. 

As discussed by Nuss (1989), convergence near the cyclone center and warm front 
ploy a significant role in the vertical circulation of the cyclone. Because of the heat flux 
distribution, stronger wind stress results to the northeast of the surface low in Case 3 
compared with Cases | and 2. The curl of the surface stress is proportional to the fric- 
tronal PBL convergence (Fleagle and Nuss, 1955), With inereased convergence into the 
surface low and warm front in Case 3 as a result of the stronger heating. the PBL 
ageostrophic flow toward the front and low is strongest in this case. This results in a 
stronger low-level transverse component of the secondary circulation, which in turn en- 
hances the intensity of the cold convevor belt flowing towards the cyclone ahead of the 
warm front, which transports boundary laver air that has been heated and moistened 
by the surface heat fluxes toward the cyclone center. The result is increased vertical 
transport of heat and moisture out of the boundary laver, contributing to stronger 
midtrepospheric heating and greater intensification of the cyclone than in Case I. In 
both Cases | and 3, the strongest ascent was found to be along a line extending from the 
cyclone center southeast and along the warm front. In Case 3, both the frictional con- 
vergence and latent heat release are offset by the frontolytical effect of the surface heat 
flux distribuuon after 24 h. It appears that the reduction in the deepening rate of the 
c\clone can be related to the removal of the favorable conditions of the first 24 hours. 
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It remaims to be shown whether the stronger frictional convergence in Case 3 is due to 


PBL processes or stronger geostrophic forcing that results from other indirect effects of 
the surface fluxes. 

The role of surface heat and moisture fluxes in oceanic cyclogenesis is complex. The 
frontogenetical effects of surface heat fluxes on the PBL stratification and frictional 
convergence are opposed by the direct frontolvtical effects on the thermal! gradient. On 
the other hand. turbulent mixing can transport heat and moisture out of the boundary 
Javer into the troposphere where it can enhance the temperature gradient through the 
release of latent heat. This feeds back to the cyclone. since the enhanced thermal gra- 
dient enhances the secondary circulation and intensity of the evclone, which in turn en- 
hances the winds and the geostrophie frictional convergence. It can be seen that the 
advance of the cyclone into a new distribution of differential heat and moisture fluxes 
has a critical nnpact on the intensity of cyclogenesis. Case 1 allows an analvsis of 
boundary laver processes in the presence of a continuously changing cross-front SST 
distribution. Case 2 allows analvsis of frontogenesis and cyclogenesis in the absence of 
surface heat arid moisture Muses. and Case 3 allows analvsis of these processes in the 
presence of an intually uniform SST gradient followed by a changing gradient after 24 


Hours ot vey cloniment. 


De TOTAL FORCING AT THE SURFACE 

The intensity of the surface warm front can be represented by the total geostrophic 
Oi acsostraphie forcing ut the surface. Pigs. 12 and 13 show the total forcing at the 
surface at 24h for Cases 1 and 2, and Fig. 14 shows the same for Case 3. Recalhng the 
direcuon of flow around verucally-ortented centers of positive and negative forcing as 
discussed in Section I]. one can determine the sense of flow in the vertical plane from 
tie sign of the forcing at the surface. In each case there is a region of shghtlv negative 
forcing ahead of the warm front, which imphes clockwise circulation in the vertical plane 
and conseyuently flow towards the northeast m the horizontal plane of the surface. 
Similarh. the regions of positive forcing imply counterclockwise circulation in the verti- 
cal plane. and consequently surface flow towards the southwest. The vertical plane of 
Sow is assumed to be perpendicular to the axis of horizontal forcing, and the sign of the 
forcing at the surfuce 1s assumed to extend Vertically above that region. 

The boundarics between the positive and negative regions of forcing are of particular 
portance since they represent the transition from counterclockwise to clockwise flow, 
and hence are regions of upward or downward vertical motion. For example, in Case 


3. the boundary between the negative region to the northeast of the front and the posi- 


ty 
"Ay 


tive region within the surface front 1s a downward vertical motion region. Similarly, the 
boundary between this positive region and the negative region south of the surface low 
is an upward vertical motion region. This combination of downward motion ahead of 
the front and upward motion behind the front supports the direct circulation and 1s 
frontogenetical in that the forcing enhances the temperature gradient. The strongest 
surface flow is to the southwest, perpendicular to the axis of forcing. which 1s close to 
the position of the warm front. 

The total forcing is strongest in Case 2 and weakest in Case 1, which supports the 
results of Nuss (1989), who found that the Case 2 warm front was the strongest. Yet this 
case did not have the strongest vertical circulation or cold conveyor belt, nor did it have 
the most rainfall, although it had the most intense evclone. The structure of the vertical 
circulation and forcing above the surface is not necessarily the same as that of the sur- 
face. A clearer picture can be obtained from an analysis of the cross-sections of total 
forcing for cach case. Of more importance is how each component of forcing is con- 
tributing to the total forcing at all levels. and what effect the surface fluxes have on these 


components. 
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Fig. 12. 
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Total forcing at the surface at 24 h, Case 1. Geopotential height in m 
(solid), total forcing in K 100 kyr! day! (dashed), and observed wind 
vectors inm s '. Wind vectors vary in length according to speed, with the 
maximum length 50 m s-'. Contour interval of forcing is 5. Heavy lines 


denote zero contours. 
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Fig. 13, Asin Fig. 12, except for Case 2. Contour interval is 10. 
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Fig. 14. As in Fig. 12, except for Case 3. Contour interval is 10. 
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Iv. RESULTS 


An ON BRS EV 

A series of surface plots of total frontogenetical forcing through the crvione domain 
were taken at six-hour intervals to analyze the evolution of total for-ing. Figs. 12 
through 14 are representative of the general pattern and relative intensities of forcing 
of the cold and warm fronts observed throughout the simulations. Although there is 
positive forcing of the warm front in each case, total forcing at the surface in Cases 1 
and 3 is strongest in the cold air sector behind the cold front. whereas in Case 2 total 
forcing Is strongest in the vicinity of the warm front. This suggests that processes unique 


to Cases 1 and 3 are dominating total forcing at the surface. 


B. TOTAL FORCING OF THE WARM FRONT 
1. Surface 

As seen in Figs. 12 through 14, the pattern of forcing of the warm front in Cases 
> and 3 is elliptical, Whereas a distinguishable ellipse of forcing does not become evident 
until 30 h im Case 1, Total forcing reaches a maximum at 30 h in Case 3, and 36 h in 
Cases | and 2. Fig. )3 1s 2 time Semes of the total forcing of the warm front, measined 
as the maximum contour value at the approximate mud point of the surface warm front. 
This is m the location of maximum forcing, and also avoids the occluded front. Total 
forcing reaches a maximum of 10 K 100 Avr! day in Case }, 40 K 100 Ayn day in 
Case 3, and 100 K 100 kn day in Case 2. Comparing the SLP deepening rates from 
Vig. ]] with the rates of growth and relative intensities of total forcing from Fig. 14 for 
each case, one can see that the strongest frontogenetical forcing and deepest cyclone are 


both in Case 2, and the weakest total forcing and slowest deepening rate are m Case 1. 
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hice 15: Evolution of total forcing of the surface warm front. 


2. Cross Sections 

Cross sections of total forcing of the warm front were taken at six-hour intervals 
for each case. oriented perpendicular to the warm front and bisecting the region of 
maximum total forcing of the warm front at the surface for that time. The cross sections 
depict more clearly the location of total forcing in the vertical plane with respect to the 
Vertical circulation field and surface warm front. By 12 hours all three cases have broad 
regions of positive total forcing in the cold air sector ahead of the warm front, and gen- 
erally negative forcing in the warm air sector behind the warm front. The counter- 
clockwise circulation induced by the positive forcing ahead of the front, and the 
clockwise circulation induced by the negative forcing behind the front, set up the direct 
circulation cell and upward vertical motion region in the warm front, as suggested in 
Section I] and as indicated in ligs. 16 through 18 by the vectors of agcostrophic vertical 
motion. Figs. 16 and 17 show the cross sections of total forcing at 24 hours for Cases 
1 and 2, and Vig. 18 shows the same for Case 3. The orientation from left to right in the 
cross-sections 1s northeast to southwest, with the cold air sector ahead of the front to 
the left and the warm air sector behind the front to the right. 

The initial model cyclone does not have an initial vertical circulation field; rather 
it evolves over the course of the simulation as frontogenetical processes occur in the 


model. The presence of a distinct updraft region, as depicted by the vectors of 


3] 


ageostrophic flow around the warm front, is evident by the 12 hour point in Case 3, with 
a less distinct updraft region in Case I at this time. Case 2 does not develop a distin- 
guishable updraft region until 18 hours. One can get a sense of the intensity of the cold 
convevor belt (CCB) from the magnitude of the ageostrophic vectors circulating around 
the direct cell in each case. The strongest CCB ts in Case 3, while the weakest is in Case 


4 


doe 


242 XX MARCH 1986 o—EQ EG) va 


1808 


(36 @ 38.0) (24 i) Wil lh 
® Se0e-02 
MAXIMUM VECTOR 





Fig. 16. 24-h cross section of total ageostrophic forcing, Case 1. Forcing is in Kk 
100 kar! day (solid), 0, in K (dashed), and ageostrophic circulation vec- 
tors. Ageostrophic circulation vectors vary in length according to speed, 
with the maximum length corresponding to 50 ms! in the horizontal and 
50 wb s-' in the vertical. Ileavy lines denote zero contours. Contour in- 
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Fig. 17. As in Fig. 16, except for Case 2. Contour interval is 5. 
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Fig. 18. Asin Fig. 16, except for Case 3. Contour imterval is 10, 


By the 18 hour point Cases | and 3 have distinguishable cllipses of positive total 
forcing extending vertically within and ahead of the updraft region, whereas in Case 2a 
sinular region is not apparent until 24 hours. Cores of positive forcing in the boundary 
layer are first distinguishable in Case 1 at 30h, and at 24 h in Cases 2 and 3, after which 
the boundary layer remains the level at which total forcing ts the strongest. Of all the 
cases, Case 2 develops the strongest boundary layer total forcing, reaching a maximum 
of 150 K 100 Ana! day at 36h. The relative maxima and times of occurrence in Cases 
Land 3 are 35 K 100 Avr! day! at 36h and 60 K 100 Ant! day! at 30 bh respectively. 
Tig. 19 shows the evolution of total forcing in the PBL for each case, measured from the 
respective cross sections as the maximum positive value in the core of positive forcing 
near the updraft region, typically at about 950 mb. Althougl: this is a crude measure 
of the frontogenesis, it gives at least qualitative evidence of how the forcing evolves m 
each case. Case 3 has the strongest initial forcing, but is the case mm which the PBL 
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forcing first declines. Between 18 bh and 24h the total forcing in Cases 2 and 3 is thre 
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same, after which the rate of growth begins to decline in Case 3 while the growth rate 


in Case 2 continues to increase. 
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Fig. 19. — Evolution of total forcing in the PBL. 


At mid levels, Case 1 develops the strongest total forcing, reaching 35 K 100 
knv! day ' at 42h. Case 2 develops a maxintum of 30 K 100 ki! day-', and Case 3 
develops a maximum of 20 K 100 faa! day! at 24h. Tig. 20 is a summary of the total 
forcing at mid levels, measured in the direct circulation cell as the maximum contour of 


positive forcing at 600 mb. 
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Fig. 20. Evolution of total forcing at mid levels. 


C. TOTAL FORCING AT 24 TOURS 

For continuity with Nuss (1989), the 24-hour time was chosen for detailed analysis. 
This time 1s representative of the most rapid phase of development. and ts just prior to 
formation of the occluded front. By this time in the evolution, the effect of the fluxes 
Will have mfluenced all the terms to some degree. In all three cases there ts a large region 
of positive forcing ahead of the warm front, and generally negative forcing behind the 
front. The gradient of total forcing at mid levels is strongest in the updraft region where 
the vertical vectors are the strongest. The positive forcing in the direct circulation cell 
dominates the negative forcing. suggesting that the direct cell is the dominant cell im the 
secondary circulation of the warm front. 

The patterns of forcing of the direct circulation cell in each case, as seen in ligs. 
16 through 18, are not symmetrical ellipses as discussed in Section I], particularly in 
Cases 2 and 3 where there is a cell of strong positive forcing in the boundary laxer be- 
neath the updraft region. Case | develops a similar region by 30 hours. There are dis- 
tinct differences in the forcing at each level in the three cases, and as discussed in the 
previous section, must be a result of the differences in the surface fluxes in ea: n case. 
As described by Nuss (1989), different surface heat and moisture flux distributions re- 
sulted in differing rates of cvclogenesis, frontogenesis, and intensity of the cold conveyor 


belt. There is certainly evidence of a relationship between total forcing of the warm 
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front, especially at the surface, and the cyclone deepening rate. The important question 
is how these flux distributions, which differed primarily near the warm front. modified 
the warm frontal temperature gradient and vertical motion in each case, as indicated by 


differences in total forcing. 


D. CONTRIBUTIONS BY EACH COMPONENT OF FORCING 
1. Overview 

The total forcing in each case is strongest in the boundary laver, with confluence 
and shear the dominant contributors to total forcing at this level. Confluent forcing 1s 
generally concentrated beneath the updraft region. with shear extending from this lo- 
cation into the direct circulation cell at low levels. Diabatic heating is generally negative 
at this level in all three cases, concentrated beneath the updraft region but extending 
ahead of the warm front as with shear. Vrictional forcing is a small component of total 
forciig of the warm front. The forcing due to diabatic heating is the dominant compo- 
nent of total forcing at mid and upper levels, and generally located within and ahead of 
the updraft region. Centers of mid-level confluence and shear are located within the di- 
rect cell ahead of the warm front. Interesting departures from these gencral results were 
found in eae case, as described below. 

peeease | 

The 24-h sensible and latent heat flux distributions for Case 1] are shown in Fig. 
21, with the positions of the front and cross-section added for reference. By this time 
the inital 1009 mb low has deepened to a 986 mb low, with strong positive total (latent 
plus sensible) heat fluxes west of the low and downward fluxes east of the low. This 
distribution prevents the formation of concentrated surface temperature gradients and 
strong temperature advection, since the upward fluxes are warming the cold air region 
and the downward fluxes are cooling the warm air region. However this effect extends 


only to the top of the boundary laver (Nuss and Anthes, 1987). 
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Fig 2). 24-h heat flux distributions for Case 1. (a) sensible and (b) latent heat 


flux distributions in W oat ?. Dashed lines indicate downward or negative 


fluxes. Fronts added for reference. From Nuss, 1989. 


a. Surface 

The frictional convergence and diabatic heating components of forcing for 
the surface are shown in Tigs. 22 and 23. Frictional forcing is negative in the region of 
upward heat fluxes northeast of the surface low, but also slightly negative in the region 
of downward heat fluxes behind the warm front. Jlowever, the frictional foreing 1s 
strongly positive in the cold air scetor behind the cold front where the upward heat flux 
is the strongest. Ihe inconsistency of the frictional forcing relative to the sign of the 
heat fluxes suggests that the momentum fluxes in the model depend on other factors that 
are not determined simply from the fluxes. This may be particularly important when the 
boundary laver is close to neutral stratification. 

The sense of the ageostrophic flow induced in the region of negative forcing 
northeast of the warm front is for surface flow to the northeast, with corresponding 
downward motion at the boundary between this region and the region of positive forcing 


closer to the front. Similarly, the induced flow in this region of positive forcing is to the 
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southwest, with downward flow along the boundary with the negative region to the 
northeast closer to the warm front. The boundary between this positive forcing region 
and the region of negative frictional forcing associated with the low and warm front ts 
upward. The combined flow would tend to support the direct circulation cell as dis- 
cussed in the previous section. The overall pattern of frictional forcing resembles the 
pattern of total heat fluxes. Over the course of the simulation the band of positive fric- 
tional forcing expands in size, although the magnitude does not substantially change. 
This suggests that as the cyclone becomes more intense and the geostrophic winds in- 
crease. downward momentum fluxes begin to influence frictional forcing more than they 


did initially, resulting in reduced negative frictional forcing. 
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Pigg 24-h surface forcing due to friction, Case 1. Forcing in K 100 kar! day! 
(dashed), with geopotential height (solid) in m and observed wind vectors 


inm s', Contour mterval of forcing ts I. 
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ices 24-h surface forcing due to diabatic heating, Case 1. Contour mterval ts 
5. 


Forcing due to diabatic heating at the surface ts frontolytical to the warm 
front, with the values of forcing larger than frictional forcing. Unlike the pattern of 
frictional forcing, the pattern of this component of forcing does not resemble the surface 
heat flux pattern, because the horizontal gradient of the heating is being considered. 
What 1s evident from Fig. 23 is that the surface heating gradient tends to reduce the 
low-level temperature gradient. ‘lie boundary between the negative forcing region along 
the warm front and the positive forcing region behind the front imphes downward mo- 
tion along the warm front. Consequently, there 1s opposition of the surface warm 
frontogenesis and secondary circulation. The diabatic heating component dominates 
any frontogenetical effect of the frictional forcing, giving a net frontolvtical effect of the 


boundary laver processes at this time. 
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b.  Cross-Sections 
The cross-sections m Tigs. 24 and 25 depict the vertical extent of these 
diabatic forcing terms. Frictional forcing extends up to 900 mb im the model, at which 
point boundary layer stresses are zero. In Fig. 24, the sense of flow between the positive 
region of frictional forcing ahead of the updraft region and the negative region beneath 


the updraft is frontogencetical, supporting the leading edge of the updraft region. 
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Fig. 24. 24-h cross-section of frictional forcing, Case I. 0, in K (dashed), 
ageostrophic vertical circulation vectors in wb s-', and ageostrophic forc- 
ing ink LOO ka! day! (solid). Contour mterval is 1. Heavy lines denote 


Zcro contours. 
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Fig. 25. = As in Fig. 24, except for diabatic heating. Contour mterval is 5. 


The boundary between the negative and positive regions of diabatic forcing beneath the 
updraft region in Fig. 25 1s frontolytical, and dominates frictional forcing as mentioned 
above. However, this frontolytical component of forcing is dominated by positive values 
of shear and confluent geostrophic deformational forcing (not shown) in the boundary 
layer ahead of the front. At mid and upper levels diabatic heating generates 
frontogenesis and supports forcing of the updraft region throughout its vertical extent. 
The transport of heat and moisture from the surface to the mid and upper levels allows 
for latent heat to be released within and behind the warm front, which enhances the 
temperature gradient in the warm air sector. Because the forcing is positive it also en- 
hances the secondary circulation. Geostrophic deformation 1s also positive at mid levels, 
although diabatic heating dominates, especially near the updraft region. The diabatic 
forcing is most susceptible to error at upper levels where the horizontal temperature 
gradient is a munimum, so a further discussion of upper-level forcing 1s precluded, and 


is not the focus of this study. 


3. Case 2 


By 24 h the cyclone in this case has deepened to 984 mb. Although it has the 
strongest warm front, it has about the same accumulated rainfall as Case | for this time, 
and also has the weakest secondary circulation and CCB. The patterns of confluent and 
shear geostrophic deformation are similar to Case I and are presented for this case for 
Pererence. 

a. Surface 

Figs. 26 and 27 depict the surface forcing due to confluent and shear 
geostrophic deformation, and is representative of the typical pattern of this forcing seen 
in Case 1 after 24h. Fig. 28 depicts the forcing due to diabatic heating at the surface. 
With the absence of surface heat and moisture fluxes in this case, frictional forcing 1s 
negligible, suggesting that surface heat fluxes play a significant role in driving the 
downward momentum flux. Although there are no surface heat and moisture fluxes in 
this case, the pattern of forcing due to diabatic heating at the surface is similar to that 
of Cuse 1. with stronger values. This suggests that either the errors in the diabatic 
heating calculation are relatively large or that there are significant differences in the 
boundary laver parameterization im the model. 

The sense of how the positive and negative regions of forcing enhance or 
retard the secondary circulation is the same as in Case ]. Again, the boundary between 
the band of positive forcing near the warm front in Fig. 28 and the broad negative region 
to the northeast of the front is a downward motion region. Simtlarly the boundary be- 
uveen tls positive region and the negative region near the surface low is upward. As 


in Case fT. this 1s frontolvtcal. 
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Fig. 26.  24-h confluent geostrophic forcing at the surface, Case 2. Contour inter- 
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As in Fig. 26, ercept for forcing due to shear. Contour interval is 5. 
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Fig. 28. = As in Fig. 26. except for diabatic heating. Contour mterval ts 5. 


b. Cross-Sections 

The strength of the vertical circulation is weaker in this case than in Case 
1, as indicated by the vertical velocity vectors of Pigs. 29 through 31. Confluence and 
shear, in Figs. 29 and 30 respectively, are strongest m the boundary layer beneath the 
updraft region. In Fig. 31, diabatic heating at mid levels, in the form of latent heat, ts 
not as strong as in Case I. This ts to be expected since the only direct source of moisture 
in the model is from the PBL. Not only rs there no {ricttonal confluence, there 1s no 
moisture transported by the CCB to surface low. In addition, the weaker secondary 
circulation in this case supports the weaker CCB found by Nuss (1989). Elowever, this 
case 1s dominated more strongly by geostrophic deformation at the surface than Case 
1, which supports the findings in previous studies that geostrophic deformation domi- 
nates frontogenesis and cyclogenesis, since this case has the strongest warm front and 


eventually the deepest cyclone. 


46 


XX MARCH 1986 


182a¢@ 
(38.@ 28.0) (28.8 16.9) 
@ 500£+82 


MAX {MUR VECTOR 





Fig. 29. 24-h cross-section of confluent geostrophic forcing, Case 2. Contour in- 


terval ts 5. 
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Fig. 30. As in Fig. 29, except for forcing due to shear. Contour imterval is 5. 
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Fig. 31. Asin Fig. 29, except for diabatic heating. Contour interval is 5. 


4. Case 3 
The cyclone in this case has developed as rapidly as the Case 2 cyclone, and al- 
though the surface front is not as intense as in Case 2 it has the most accumulated 
rainfall at 24h. The heat flux distribution at 24h, as shown in Tig. 32, is characterized 
by strong upward heat fluxes behind the cold front, positive heat fluxes northeast of the 
surface low and warm front, and a pocket of downward heat Nuxes near the trailing cdge 
of the center of the warm front, similar to Case 1. 
a. Surface 
The pattern of frictional forcing in Case 3, as shown in Jig. 33, Is quite 
different to that of Case I, in that alinost the entire region east of the surface low Is 
negative. Within this negative region is a pocket of slightly positive frictional forcing 
near the warm front. Again, the frictional forcing distribution 1s more complex than the 
heat flux pattern, indicating the significance of other factors. Whereas in Case | the 


region of positive forcing was ahead of the warm front, in this case itis located behind 
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the warm front. The resultant forcing of the flow between the positive and negative re- 
gions does not seem to located in the same position relative to the warm front to en- 
hance the direct circulation cell and frontogenesis. As in Case I, this region of positive 
frictional forcing expands in size, although not in magnitude. ‘The magnitude of negative 
frictional forcing does increase slightly, although it decreases in area. This latter result 
was not seen in Case 1, and suggests that since the two cyclones are of equal intensity 
at this time, the difference is due to the surface heat flux distribution. 

The 24-h surface diabatic heating in this case, as seen in Fig. 34, is similar 
to the other cases, but with a more narrow band of negative forcing associated with the 
warm front. The strength of the gradient from the negative to positive regions is 
strongest in this case, supporting a stronger transverse flow and vertical flow at the 
boundaries. Just as frictional forcing appears to be Iess frontogcnetical than in Case 1, 
the locations of the upward and downward motion regions appear to be less frontolytcal 


than in the other two cases. An analysis of the cross sections yields a clearer picture of 


this relationship. 





Fig. 32. Asin Fig. 21, except for Case 3. 
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Fig. 33. 24-h surface forcing due to friction, Case 3. Contour imterval is 2. 
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Fig. 34. 9 24-h surface forcing due to diabatic heating, Case 3. Contour interval is 
2. 


b.  Cross-Sections 

The cross-section of frictional forcing in Tig. 35 reveals that the region of 
positive frictional forcing is located beneath the trailing edge (warm sector side) of the 
direct circulation cell updraft region, with negative forcing directly beneath and ahead 
of the updraft region. The resultant vertical motion between these centers of foreing acts 
to oppose the vertical component of the direct circulation cell. Tfowever, because the 
negative frictional forcing ahead of the front is strongest in this case compared with Case 
1, one can imagine a stronger transverse flow towards the front and cyclone at the top 
of the PBL. Thus frictional forcing is not enhancing the updraft region but is enhancing 
the transverse component of the direct circulation cel! at low levels. 

This case has stronger values of mid level latent heat than the other two 
cases at this time, as seen m Tig. 36, and indeed, this case has the most rainfall accu- 


mulation at 24h. This suggests that the region of strong upward moisture flux northeast 


of the surface front has influenced the latent heat of the mid troposphere, transported 


to the front and cyclone by the influence of frictional forcing. 
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Tig. 35. = 24-h cross-section of frictional forcing, Case 3. Contour interval is I. 
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Fig. 36. As in Fig. 35, except for diabatic heating. Contour mterval is 10. 


E. TIME SERIES OF FORCING 

The results at 24 h are a snapshot at one time, and the feedback of the surface 
diabatic effects on geostrophic forcing has not been separated. A time series of the 
contribution to total forcing by each component at low, mid, and upper levels attempts 
to separate the mfluence of each forcing term as they evolve throughout the simulation. 
Forcing Was measured as the maximum value of each term along a horizontal line across 
the direct circulation cell m the respective cross-sections of each case. The most am- 
brguous term at low levels is the frictional forcing term. Frictional forcing was measured 
as the maxmmum positive value near the leading edge of the updraft region of the direct 
circulation cell, stnce this would enhance the updraft region m this location, ‘I his de- 
piction does not take tnto account the influence of the negative frictional forcing ahead 
of the warm front. Theoretically, the clockwise circulation around the negative forcing 
center is acting to enhance the low-level transverse component of the circulation cell at 


the top of the PBL where the circulation around the region is towards the front and 
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cyclone. Forcing due to diabatic heating was measured as the maximum negative value 
at the base of the updraft region, since it was generally frontolytical in each case at this 
level and location. Confluence and shear were measured as the maximum positive values 
in the direct circulation cell within and ahead of the front. The 950 mb level was chosen 
since it is the level of maximum values for most of the terms at each tme. At mid levels, 
the G00 mb level was chosen to represent the forcing due to confluence, shear, and latent 
heat release. Mid-level forcing was found to be a maximum at this level. with confluence 
and shear generally centered within the direct circulation cell. and latent heating located 
within the updraft region. Because the upper-level diabatic forcing 1s most susceptible 
to error, the maxima were located at 300 mb, where the center of the upper-level trans- 
verse component is located in each case. For reference. the SLP deepening rate is also 
provided with the 950 mb time series in each case. 
1. Boundary Layer 

The time series of forcing m the boundary laver for Cases 1 and 2 is depicted in 
Tig. 37. and for Case 3 mn Fig. 38. In all three cases. confluent and shear geostrophic 
deformational forcing dominate total forcing at low levels, with confluence dominating 
shear throughout the simulations except for Case 1. where shear dominates confluence. 
Total forcing does not peak until 30 h in Case 3 and 36h in Cases 1 and 2. This 1s 
consistent with the decline in the deepening rate of the cyclone in each case at these 
times. Surface diabatic heating is frontolytical, initially peaking at the same time as 
evostrophic forcing. This component of low level forcing is frontolytical in all three 
cases, and reduces total forcing by about 20 °«o. Unlike geostrophic deformational 
forcing, diabatic heating at the surface does not have a single peak and period of decline. 
Frictional forcing 1s not a significant contributor to total forcing in Cases | and 3, but 
it is the first component of total forcing to develop in Case 3. In all three cases, 
confluent and shear geostrophic deformational forcing and diabatic heating do not de- 
Vela ceelow lewels until alter Wh in Case 3 and 18 h in Cases 1 ande2. Thiec earlier de- 
velopment of geostrophic deformation in Case 3 is presumably due to the influence of 
the frictional forcing that was present from the beginning. The period of maximum de- 
velopment of geostrophic and diabatic forcing is between 1§ h and 24 h in all three cases, 


and the period of decline 1s between 30 h and 36h. 
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Fig. 37. ‘Time series of forcing in the PBL. (a) Case I, and (b) Case 2. Note the 
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Fig. 38. —As in Fig. 37, except for Case 3. 


2. Mid and Upper Levels 

The time series of total forcing at mid levels 1s depicted in Figs. 39 and 40, The 
time scries for upper Ievels in depicted in Pigs. 41 and 42. At both levels, diabatic 
heatmg dominates total forcing. except for Case 3 where geostrophic deformational 
forcing dominates. Mid level diabatic heating in the form of latent heat release ts 
strongest in Case 1, where it docs not peak and level off until 42 h. twelve hours after 
the deepening rate has begun to decline. In Cases 2 and 3 dinbatic heating reaches a 
steady state at about 24 hous. At upper levels, there ts some contribution to total 
forcing by confluent and shear gcostropluc deformation. In Cases | and 2, diabatic 
heating at this level peaks at 36 h. then declines, whereas in Case 3 it remains steady 


alter 24 h. 
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Fig. 39. = Asin Fig. 37, except for 600 mb. 
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Fig. 40. Asin Fig. 38, except for 600 mb. 
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Fig. 42. As in Fig. 38, except for 300 mb. 


F. COMPARISON OF CASES 
1. Boundary Layer 
Case 2 is similar to Case | during the first 18 hours, after which time the most 
significant differences are the rapid rise in geostrophic forcing in Case 2, the dominance 
of shear over confluence in Case 1, and the larger surface frontolytical effects of diabatic 
heating in Case 2. Comparing Case 3 with Case I, the most significant differences are 
the stronger frictional forcing in the first 12 hours of development, the rapid rise in 
gcostrophic forcing similar to Case 2 for the first 24 hours, and a decline in the 
frontolytical effects of diabatic heating after 30 hours. 
2. Mid and Upper Levels 
The most significant differences between Cases | and 2 are the steady rise in 
diabatic heating in Case | coupled with the presence of geostrophic forcing. ‘There is 
no significant geostrophic forcing at mid levels in Case 2. In Case 3, latent heating is 
evident, but it not the dominant mid level forcing as it is in Case 1. Both confluent and 
shear deformational forcing dominate diabatic heating in Case 3. At upper levels, the 
only significant difference is the absence of any significant geostrophic deformational 


forcing in Case 2. 
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V. DISCUSSION 


As OVERVIEW 

The basic conceptual model of frontogenesis is that geostrophic deformation forces 
strong frontogenesis at surface and upper levels, and at mid levels frontogenesis is usu- 
ally weak and forced primarily by latent heat release. Frontogenetical forcing can be 
modified by small-scale frictional forcing and diabatic heating, with their relative con- 
tributions to frontogenesis modified by surface and boundary layer processes. 

A review of the SLP deepening rates from Fig. 11 shows that for the first 12 hours 
each case deepens at the same rate, indicating that the surface flux distribution peculiar 
to each case does not have an apparent effect on frontogenesis until 12 h. None of the 
cases has significant geostrophic deformation enhancing the forcing during the first 12 
hours, suggesting that small-scale processes of frictional forcing and diubatic heating in 
the boundary laver are important in setting the stage for the rate of development and 
structure of the cvclone and warm front. Once cvclogenesis has begun, forcing due to 
geostrophie deformation builds, which then feeds back on itself and dominates the fore- 
ing. In later stages, friction and diabatic heating act to modify the geostrophic defor- 


matronal forcing. 


B. BOUNDARY LAYER 
1. Geostrophic Deformational Forcing 

In their two-dimensional model simulation of frontogenesis. Baldwin et al. 
(1984) found that geostrophic deformation accounted for the large-scale features of the 
secondary circulation around a cold front, but was not the dominant forcing at low lev- 
els. Bond and Hleagle (1985) found similar results in their case study of an actual 
oceanic cold front. Baldwin et al. (1984) found that shear deformation dominated con- 
fluence in their model simulation. With respect to the forcing of a warm front, this study 
found that geostrophic deformation was donunant at low levels, with confluence domi- 
nating shear im those cases with the most intense warm fronts. Mudrick (1974) found 
that confluent geostrophic deformation was favorable for warm frontogenesis. with 
shear geostrophic deformation favorable for cold frontogenesis. 

‘As mentioned in the last section, confluence dominates shear in Cases 2 and 3, 
and shear geostrophic deformational forcing dominates confluence at all levels in Case 


1. Shear deformation is not as efficient for frontogenesis as confluence aud may be a 


partial explanation for the lower intensity of the warm front and cyclone in Case 1, The 
strong frontolvtical effect of the surface fluxes in Case 1 to reduce the temperature gra- 
dient also prevents the development of a frontal pressure trough. Lack of significant 
curvature precludes the development of strong confluence. This contrasts with Cases 2 
and 3. which develop stronger temperature gradients and more of a pressure trough 
along the front. In Case 3, which is the first case to develop significant gesstrophic de- 
formational forcing. the frontogenctical effect of the friction apparenth helps to con- 
centrate the temperature gradient more rapidly in this case. This occurs in spite of 
relatively strong frontolvsis due to surface heating. 
2. Frictional Forcing 

Baldwin et al. (1984) found that frictional forcing was the dominant forcing at 
low levels, and Bond and Fleagle (1985) found that frictional convergence in the 
boundary laver accounted for 90°o of the observed total forcing at this level. The results 
of this study suggest that frictional forcing is not nearly as strong a component of forc- 
ing of the warm front, accounting for less than five percent of forcing at low levels. In 
addition. the pattern of frictional forcing ahead of the warm front is not as consistent in 
its contribution to foreing of the warm front as it 1s to the forcing of the cold front found 
in these other studies. It is not known whether this weaker forcing by friction in the 
warm front studied here is due to model characteristics or some fundamental difference 
between warm and cold fronts. The stronger frictional forcing and more intense cold 
front in Case 3 compared with Case 1 suggest that frictional forcing can greatly enhance 
frontegenesis. This enhancement is partially due to indirect effects on other more 
dominant processes. 

In addition, frictional forcing of the warm front 1s not consistently predictable 
as to whether it is frontogenctical. As shown im the results. sometimes the boundary 
between the regions of positive and negative forcing was such as to oppose the forcing 
by the other terms and sometimes it seems to enhance the foremg by the other terms. 
The development of these boundaries seems to be closely related to the position of the 
cyclone and warm front with respect to the surface fluxes, but apparently depends on 
other factors as well. Surface fluxes remain strongly positive behind the cold front, 
Whereas the flux distribution across the warm front 1s confused by the pocket of down- 
ward heat fluxes near the front and the track of the cyclone relative to the SST gradient. 

There is the suggestion of a trme lag between frictional forcing and its effect on 
geostrophic deformational forcing. in particular confluent forcing, in Case 3. The initial 


development of confluent forcing is preceded by 12 hours of steady, positive frictional 
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forcing ahead of the warm front, and geostrophic forcing develops 6 hours prior to 
similar forcing in the other two cases, which have little or no frictional forcing. Simi- 
larly, the decrease in the rate of growth of the confluent (and shear) forcing is preceded 
by 12 hours of increasingly negative frictional forcing in this region. The reversal in the 
frictional forcing mav be influenced by the track of the cyclone relative to the SST gra- 
dient. During the first 24 hours, the cvclone tracks towards the region of strong positive 
heat fluxes and the frictional forcing is frontogenetical. When the cxclone begins to 
move into colder water and away from the influence of the positive heat fluxes, the 
frictional forcing becomes frontolytical. 
3. Diabatie Heating 

Baldwin et al. (1984), using a thermally insulated boundary laver in their model, 
found that boundary laver heat fluxes had a minimal effect on the svnoptic-scale sec- 
ondarv eirculation. Bond and Fleagle (1985) included boundary laver heat fluxes in their 
study, yet found similar results, although boundary layer turbulent heat fluxes were 
found to be frontolytical in the boundary laver. Nuss and Anthes (1987) found that 
boundary laver heat fluxes influenced the structure only below 850 mb, suggesting that 
surface heat flux effects are confined to the boundary layer, due to the imflucnce of ihe 
mversion laver above the PBL. This study found similar results in that forcing due to 
diabatic heating 1s frontolvtical and occurs at low levels only. Its magnitude is compa- 
rable to that of shear deformutional foremg in Case 2 and confluent forcing in Case 1, 
and in Case 3 itis about half that of shear. Of interest 1s the presence of negative forcing 
in the PBL in Case 2 despite the absence of surface heat and moisture fluxes. This sug- 
gests that in the absence of differential surface fluxes. there is surface cooling perceived 
by the model. This may be the wav the model handles convective versus non-convective 
precipitation, or it may be a result of the computation of this term of forcing due to 


advection. 


C. MID AND UPPER LEVELS 
1. Geostrophic Deformation 
As with studies of cold front forcing. geostrophic deformational forcing of the 
warm front ut nud levels is dominant within the direct circulation cell ahead of the front, 
and at upper levels it enhances the upper-level transverse component of the secondary 
circulation. THflowever, the overall values are typically about 25°o of the low-level 
geostrophic forcing. The weakest mid- and upper-level geostrophic forcing is in Case 2, 


suggesting that mid- and upper-level geostrophic forcing is not as criucal to the devel- 
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opment of the cyclone as it is at low levels, since this case has the deepest cyclone and 
Strongest geostrophic forcing at the surface. 
2. Diabatic Heating 

Diabatic heating in the form of latent heat release was found to be the dominant 
forcing mechanism for the vertical component of the secondary circulation, with positive 
centers located within and slightly ahead of the updraft region and warm front. All three 
cases develop some mid-level diabatic forcing due to latent heating in ‘ne model. The 
strongest nud-level latent heating developed in Case | and it continued to merease until 
42h. This case was the weakest cvclone. which indicates that latent heating and a strong 
vertical circulation in the warm front are not the only factors for strong cyclogenesis. 
Case 3 developes a relatively steady-state of latent heating bv 24 h. and the strongest 
vertical motion, however the diabatic forcing was less than Case I. This is presumably 
a feature of the cumulus parameterization in the model. Case 2 has the deepest cyclone 
by 24h but produced the weakest mid-level diabatic forcing. This lack of latent heating 
is probably a reflection of the lack of moisture input from the surface in this case. These 
results suggest that details of the cumulus parameterization and moist processes in the 
model strongly alter the intensity of the mid-level diabatic forcing of frontogenesis and 
the secondary circulation. Since the diabatic forcing is dominant at mid and upper lev- 
els, the circulation at these levels will be sensitive to model errors in precipitation proc- 


esses, 


VI. CONCLUSIONS 


The results of the model simulations confirm the results of previous observational 
and numerical studies of frontogenests. They also suggest important interactions be- 
tween the surface heat and moisture flux distributions, the forcing of the warm front 
secondary circulation, and cvclogenesis. Geostrophic deformational forcing was domi- 
nant. but evolved differently depending upon surface friction and diabatic heating influ- 
ences, 

The presence of frictional forcing enhances the low-level transverse component of 
the secondary circulation, which in turn enhances the cold conveyor belt. With strong 
positive heat fluxes in a favorable location, as in Case 3, frictional forcing appears to 
enhance development of geostrophic deformational forcing at an earlier stage of 
cvclogenesis, which also enhances the secondarv circulation. As the cvclone and 
geostrophic deformational forcing intensify, the surface forcing due to friction acts only 
to modify geostrophic forcmg. The development of positive frictional forcing in the 
latter stages of cvclogenesis, as in Case |. had hittle effect on the decpenine rate, except 
perhaps to maintain the geostrophic forcing at such a level as to vield the final 4 mb 
drop before cyclogenesis finally ceased. 

The effect of surface heat and moisture fluxes and their distribution on forcing of the 
secondary circulation and cyclogenesis is more complex. and the results less definitive 
than with frictional forcing. Surface forcing due to diabatic heating is consistently 
frontolvtical in each case. The fact that it 1s a strong component of total forcing in the 
absence of surface heat and moistrure fluxes, as in Case 2, suggests that either the error 
ts large or there are significant boundary laver differences in the model. Early surface 
relaxauon of the frontal temperature gradient as in Case | appears to reduce the poten- 


tial deepening rate of the cyclone, and dominates frictional forcing. 
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VII. RECOMMENDATIONS 


The results suggest areas for further research, including: 


Verify the magnitudes of the various forcing terms in an actual case study. Because 
the model is drier than real oceanic cyclones, forcing due to diabatic heating effects may 
be different. 

Verifv’ the effect of surface fluxes and boundary laver structure on the frictional 
forcing of the warm front. 

Better understand the role of various model parameterizations involved in deter- 
nuning convective versus non-convective release of latent heat. This may help to explain 
the presence of negative surface diabatic forcing in the absence of surface heat and 


moisture fluxes. 
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